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We compute the O(a^) contributions to the photon energy spectrum of the inclu- 
sive decay B Xg'^ associated with the magnetic penguin operator O7. They are 
an essential part of the ongoing NNLO calculation of this important decay. We use 
two different methods to evaluate the master integrals, one based on the differential 
equation approach and the other on sector decomposition, leading to identical results 
which in turn agree with those of a recent independent calculation by Melnikov and 
Mitov. We study the numerical relevance of this set of NNLO contributions in the 
photon energy spectrum and discuss the change of bottom quark mass scheme. 
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Abstract 



1 Introduction 



More than a decade after their first direct observation, radiative B decays have become a key 
element in the program of precision tests of the Standard Model (SM) and its extensions. 
The inclusive decay B ^ Xg^ is particularly well suited to this precision program thanks to 
the low sensitivity to non-perturbative effects. The present experimental world average [1] 
for the branching ratio of S — > 7 has a total error of about 6% and agrees well with the 
SM prediction, that is subject to a considerably larger uncertainty [2]. In view of the final 
accuracy expected at the B factories, about 5%, the SM calculation needs to be improved. It 
presently includes ncxt-to-lcading order (NLO) perturbative QCD corrections as well as the 
leading non-perturbative and electroweak effects (see e.g. [3] for a complete list of references) . 
The calculation of next-to-next-to-leading order (NNLO) QCD effects is currently under way 
and is expected to bring the theoretical accuracy to the required level. 

Among the NNLO QCD contributions, the following have already been computed: i) all 

relevant 0{a1) Wilson coefficients [4]; ii) the relevant three loop anomalous dimension matrix 
[5,6]; iii) the 0{alNF) corrections to the matrix elements of the operators Oi, 02,07,08 [7]; 
iv) those matrix elements related to the dominant electromagnetic dipole operator O7 [8-10]. 
Among the pieces still missing to date, there are two rather important and difficult ones. 
The first concerns the four-loop anomalous dimension matrix elements that describe the 
0(ag) mixing of the four-quark operators into O7 and Og; preliminary results on this have 
been recently presented [11]. The second concerns the three- loop, 0(q;^), matrix elements 
of the four-quark operators Oi and O2. These loops contain the charm quark, and the 
ambiguity associated to the choice of scale for its mass is the main source of uncertainty at 
NLO [2]. A NNLO calculation of these matrix elements is therefore crucial to reduce the 
overall theoretical error on the branching ratio [12]. Other missing contributions, like the 
0{af) matrix elements of Og, may also prove essential to reach a high theoretical precision. 

In addition to the total branching fraction, the photon energy spectrum in 5 ^ 7 is 
also a useful observable: it receives both perturbative and non-perturbative contributions 
and is important for several reasons: i) the precise knowledge of the spectrum is necessary to 
predict the measured branching ratio since experiments apply a lower cut of 1.8-2.0 GeV on 
the photon energy in the B center-of-mass frame; ii) it is almost insensitive to new physics; 
iii) from the moments of the truncated spectrum one can extract relevant information on 
the parameters of the Heavy Quark Expansion. In particular, a precise value of the h 
quark mass can be extracted from the mean value {E^) of the photon energy [13, 14]; iv) 
the measured spectrum gives direct information on the Shape Function that encodes the 
QCD dynamics in the endpoint region, where non-perturbative contributions dominate and 
perturbative corrections must be resummed [15]. The NNLO Sudakov resummation was first 
completed in [16] and there has been recent progress towards NNLO in the multi-scale OPE 
approach [17]. A detailed knowledge of the Shape Function is useful for the determination 
of \Vub\ from inclusive semileptonic B decays. 

The perturbative contributions to the photon energy spectrum are known at 0{as) [18] 
and 0(q;^/3o) [19] since several years. In this paper we present a calculation of the 0(q;^) 
photon spectrum induced by the magnetic operator O7. This operator gives the dominant 
contribution to the spectrum at 0{as)- Our results, which were obtained by using two 
different methods, confirm those of a recent paper [10]. 
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Our paper is organized in the following way. In Sec. 2 we present our analytical results for 
the unnormalized photon energy spectrum and compare with the literature. In this section 
we also compare numerically the NLO and the NNLO spectra, and study the impact of a 
change of the bottom quark mass scheme. Sec. 3 is devoted to a detailed description of the 
techniques used for our calculation. Finally, we give a brief summary in Sec. 4. 



2 Results and applications 

2.1 Analytical results for the photon energy spectrum 



Within the low-energy effective theory the partonic b — > Xg'j decay rate can be written as 

327r4 



rib - Xr^-^U^,,, = Gla^^mMml ^ ^^^^^^^ ^^.l) 

where rrib and mf,(/x) denote the pole and the running MS mass of the b quark, respectively, 
Ci^ilJ') the effective Wilson coefficients at the low energy scale, and Eq the energy cut in the 
photon spectrum. Here we give the result for the function G^riEo, /i) corresponding to the 
self-interference of the electromagnetic dipole operator 

O7 = mbii^) (sa'^^PRb) F^, (2.2) 

including 0{a^) terms as required for NNLO accuracy. Introducing the dimensionless vari- 
able 

2-E'^ 

z^—^, 2.3 

this function can be written as an integral over the (rescaled) photon energy spectrum 
dGr7{z, ij)/dz, i.e., 

Grr{Eo,f^) = f^^^^dz, (2.4) 

where z = Zq corresponds to E^ = Eq. When working to the required precision, the photon 
energy spectrum can be separated into three different parts, 

dGr7{z,n) _ ^ dG]^''{z,n) 

d; -z- Tz ' ^^-^^ 

n=2 

corresponding to the n particles in the final state, namely the fo — > 57 (n = 2), 6 — > sjg 
(n = 3), 6 — > s^gg and b s^qq (both n — A, q & {u,d,c,s}) transitions. Furthermore, 
each individual contribution can itself be written in the form 

^^Z^^ =/,(/.) 5(1-.), 
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Consequently, the complete result for the photon energy spectrum is of the form 
dGr7{z,iJ,) 



dz 



F{n) 5(1 -z)+ R^{z, n) + R^{z, n) , 



(2.6) 



(2.7) 



with -F(/u) = /2(At) + /3(Ai) + /4(Ai)- We stress here that the calculation of dG7r{z, n)/dz boils 
down to the determination of R^{z,ii) and R4{z,iJ,), because the coefficient F{iJ,) in front of 
the delta-function can then be fixed by the requirement that the total integral 
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yields the result obtained in [8] , which was recently confirmed in [9] . 

As we discuss the details of the calculation in Sec. 3, we immediately present our final 
result for the photon energy spectrum (including also the order q;° and a\ pieces): 
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where 
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The numerical values of the color factors are Ci? = 4/3, = 3, and Tr = 1/2. Furthermore, 
Nl and Nh denote the number of light (m^ = 0) and heavy (m^ = mf,) quark fiavors, that 
is the total number of quark fiavors is Np = + Nh- The functions H'^^^ and H^"^'^^ 
{j — a, na, NL, NH) appearing in (2.9) and (2.10) are given by 
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ln(/i/mfe), IjI^^z) = dxLi2{x)/x, ({3) is the Riemann zeta-function, and 
[ln"(l — x)/{l — a;)]+, with n = 0, . . . , 3, are plus-distributions defined in the standard way. 

We also worked out the normalized photon energy spectrum 1/Gn{0, fj,) ■ dG77{z, fj,)/dz, 
and after setting /i — rrib in our result, we find complete agreement with a recent paper by 
Melnikov and Mitov [10]. 



2.2 Comparison of NLO and NNLO results 

The left frame of Fig. 1 shows the dependence of the spectrum dG{z, mh)/dz on the rcscalcd 
photon energy z at NLO (dashed curve) and NNLO (solid curve). Also plotted is the NLO 
contribution supplemented by those NNLO terms which are proportional to a'^Po (dotted 
curve). These terms, already worked out in [20], are often called Brodsky-Lepage-Mackenzie 
(BLM) terms [21]. They arise from the contributions ~ a'^Ni, after replacing Nl —3f3o/2 
according to the procedure of naive non-abelianization [22]. As seen in the figure, the BLM 
terms provide the dominant part of the 0{a'l) corrections. As pointed out in [8], one should 
stress that this statement refers to the scheme where the decay width is written as in (2.1), 
i.e. with the combination m^(//)m^ explicitly appearing. 

The behavior of the non-BLM corrections was discussed in detail in [10] and since we 
agree with their findings, we do not repeat this discussion. However, when comparing the 
left frame of our Fig. 1 with Fig. 1 of [10], the reader should bear in mind that the spectrum 
shown in the latter is normalized to the fully integrated quantity 6^77(0, m^), followed by a 
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Figure 1: The spectrum dG-n^z^rrii^) / dz as a function of z (left) and ^^^jdGjj / dz\dz as a 
function of Zq (right) at NLO (red dashed curve) and NNLO (blue solid curve) for = 4, 
Nh = 1 and as{mb) = 0.22. The black dotted curve is the BLM approximation of the NNLO 
result. 



consistent expansion in powers of a^- As a result the NLO and BLM curves are equal in 
both figures, but the NNLO approximation lies below the BLM curve in the left frame of our 
Fig. 1, while it lies above the BLM approximation in the case of the normalized spectrum 
shown in Fig. 1 of [10]. 

In the right frame of Fig. 1 we display the quantity J^^ldGn / dz]dz as a function of the 
(rescaled) photon energy cut-off Zq. As can be seen, the difference between the NNLO and 
the BLM result is larger than in the left plot. This is due to the fact that the endpoint 
contributions at 2; = 1 enter this integrated quantity, for which the difference between the 
BLM approximation and the full NNLO contributions is somewhat larger. But still the BLM 
terms provide a good approximation also in this case. 

2.3 Change of bottom mass scheme 

The above results have been obtained for a pole b quark mass, that is known to be affected 
by a leading infrared renormalon and generally leads to slowly converging perturbative ex- 
pansions. In order to consider a change of scheme for this parameter, it is convenient to use 
the normalized photon energy spectrum, namely 

/(., f,) ^ = S{1 -z) + ^C,f('^ (z) + (^) ' /.) + 0{al) (2.17) 

G'77(0,/i) vr V y 

A redefinition of the b quark mass mf — Sm^ , where is the b mass in the scheme 

X, leads to a rescaling of the variable z z^{l + Sm^ /m^). A meaningful comparison 
of different mass schemes can be made at the level of moments of the spectrum or fraction 
of events. They are physical observables and are expected to have small and computable 
non-perturbative corrections. Several calculations of these quantities exist [14,23,24], but 
they do not always include all the information coming from the 0(q;^) calculation of the 
spectrum. In the following we drop the fi dependence from all formulas, set fi = m^, and 
use as{mf,) = 0.22 in the numerics. 
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Let us consider for instance the fraction R{Eq) of events with photon energy above a 
cut Eq. Neglecting as usual contributions from operators other than O7, the perturbative 
expansion for this quantity, 



R{Eo 



1 — / dz f{z), where Zq 
Jo 



2Eo 



(2.18) 



depends on the mass scheme adopted. In particular, the integral above receives contributions 
that start at 0{as). In the pole mass scheme with = 4.8 GeV and Eq — 1.8 GeV one 
obtains 



i?(1.8GeV)|poie 



0.0145 - 0.0085 -0.0019 + 



0.9751 + 



(2.19) 



where we have listed separately the trcc-level component, the NLO term, the BLM result, 
and the non-BLM terms at NNLO. As typical in b decays, the large value of the BLM term 
suggests an optimal scale for ag of 0(1-2 GeV). The inherently large uncertainty in the 
determination of the pole mass implies a significant uncertainty in R{Eo). For instance, 
using rub = 5.0 GeV we obtain i?(1.8GeV)|poie = 0.9811. 

Let us repeat the same calculation in a different scheme X for the b mass. Apart from 
the change in the numerical value of the mass in the evaluation of Zq, there is an extra NNLO 
(non-BLM) term 



X 



SE'^iEo)^ -Cfz^ 



5m. 



m, 



xf''\^^) 



(2.20) 



To illustrate the numerical difference due to a change of scheme in R{Eo), we consider here 
the kinetic mass of the bottom quark [25], also employed in [14]. At 0{al) the kinetic mass 
is related to the pole mass by [26] 



o^sirrib) ^ /4 



TT 



^"V3 2 rUb 



2 f^ku 



2mft 



(2.21) 



Here /Xkin is a Wilsonian cutoff that factorizcs soft and hard gluons; the standard choice is 
//kin = 1 GeV. Inserting (2.21) in (2.20) and employing m5"^(l GeV) = 4.6 GeV, as suggested 
by recent fits to the semileptonic moments [13], we obtain 



R{1.S GeV)|kin = 1 - 0.0196 - 0.0118 - 0.0025 + 0.0042 + • • • = 0.9703 + 



(2.22) 



where the last contribution has a positive sign and comes from (2.20). The result in the 
kinetic scheme is compatible with the one in the pole scheme (2.19), that is expected however 
to have larger higher order corrections. Eq. (2.22) can be compared with the simple estimate 
given in [10]: the non-BLM correction has the same magnitude as quoted in [10], but differs 
for the sign once the additional shift (2.20) is taken into account. As suggested in [10], one 
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can write the perturbative expansion in terms of at a scale lower than m^, improving 
the apparent convergence of the series, and increasing (by up to a factor ~ 2.5) the relative 
importance of the non-BLM term. A consistent implementation of the kinetic scheme has 
been presented, up to 0{a'^f3o), in [14]: besides the use of the kinetic bottom mass, it involves 
a proper definition of the non-perturbativc matrix elements of higher dimensional operators 
and a cut on soft-gluon radiation that modifies the shape of the spectrum close to the end- 
point. A complete implementation of these effects is beyond the scope of the present paper, 
but we expect (2.22) to show the likely impact of the 0{a^) photon spectrum on the kinetic 
scheme calculation of i?(£'o)■ 
Finally, we consider the truncated first moment of the normalized spectrum in the pole 
and kinetic mass schemes, neglecting all non-perturbative corrections and the contributions 
from operators other than O7. It can be written in the on-shell scheme as 



1 - I'dz (1 - z) fiz) (^1 + f{y)^ + 0(af ) 



(2.23) 



Using the same numerical inputs as above, the pole scheme expansion is 

{Ej)e-,>i.8 GeV 0.0246 - 0.0242 + 0.0037 - 0.0004] = 2.291 GeV, (2.24) 

where we have listed separately the tree-level, NLO, BLM, NNLO non-BLM, and the contri- 
bution of the interference (last) term in (2.23). The relative weight of the BLM contribution 
is even larger than in the previous example, as expected since the first moment probes a 
region of smaller gluon energies. 

A change of bottom mass scheme in the first moment leads to the extra term 




4)f^'\4)+ fdz{l-z)f'\, 

•J Zn 



(2.25) 



l-Cp—[z^{l-z^)f^'nz^)+ I dzil-z)f^'>{z) 

where 5m^ now includes the O(a^) term, and the r.h.s. must be expanded up to 0(al). 
The scheme change affects already the NLO calculation and at NNLO contributes also to 
the BLM term. The numerical result in the kinetic scheme is 

'2 

+0.0293 + 0.0296 - 0.0035 - 0.0016] = 2.324 GeV, (2.26) 



(^7)£^>i.8 GeV - — ^-^ ^[1 - 0.0234- 0.0234 + 0.0038-0.0005 



where the last four entries are the NLO, BLM, NNLO non-BLM mass shifts, and the mixed 
term in (2.25). Once again, the result is consistent with the one in the on-shell scheme, 
that is however subject to a larger uncertainty; it is also consistent with the experimental 
value by the Belle collaboration [27]: {E^)e^>i.s GeV = 2.292 ± 0.026 ± 0.034 GeV. We stress 
that our calculation in the kinetic scheme is not complete and that non-perturbative effects 
have not been included (see [14]). We agree with [10] for the part included in that paper, 
which does not properly consider the change of scheme. Our result for the non-BLM NNLO 
contribution to (2.26) is mf''{l GeV)/2 (0.0038-0.0005-0.0035-0.0016) ~ -4 MeV. Even 
if we rescale it to account for a lower scale, this has the opposite sign and is much smaller 
than the estimate in [10]. It can be expected to shift the h quark mass extracted from the 
first moment by less than 15 MeV, well beyond present sensitivity. On the other hand, it 
is clear that a complete NNLO analysis of the i? ^ X, 7 moments based on the spectrum 
presented here will have a smaller perturbative uncertainty. 
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Figure 2: Diagrams contributing to Ti^n for n = 3 and 4 at O(a^). We draw here only the 
nonabelian diagrams. Thick lines represent the massive b quark, thin lines the massless s 
quark, wavy lines photons and curly lines gluons. Dashed lines separate the original Feynman 
diagrams which enter the squared amplitude |M„p. The dashed blue lines indicate cuts 
contributing to the b — > S'ygg process, and the dashed orange lines indicate cuts contributing 
to the b s'-fg process. Possible left-right reflected diagrams are not shown. See text for 
more details. 

3 Details of the calculation 

Our goal is to calculate the NNLO corrections to the quantity dG-jji^z, jj) / dz introduced in 
(2.4), and, as stated in the previous section, to do so we only have to calculate contributions 
with three and four particles in the final state. Hence, we restrict the following discussion 
to the cases n = 3 and 4 [see (2.5)]. 

We start by considering the general expression for the decay rate of a massive b quark 
with momentum pi, into tt, = 3,4 massless final-state particles with momenta /cj. 




where the squared Feynman amplitude |M„p is always understood as summed over final 
spin-, polarization- and color states, and averaged over the spin directions and colors of the 
decaying b quark. For n = 4 it also includes a possible factor of 1/2 if there are two identical 
particles in the final state. Furthermore, d = A — 2e denotes the space-time dimension that 
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Figure 3: Same as in Fig. 2 for diagrams with a closed light quark loop. Similar diagrams 
with a closed gluon loop are also present. The dashed blue lines indicate cuts contributing 
to the b S'-yqq (q G {u, d, a, s} ) and b — >• sjgg processes. Concerning diagrams with 
closed ghost loops, see text. Note that three-particles cuts of the diagrams above are zero 
in dimensional regularization. 




Figure 4: Same as in Fig. 2 for diagrams including a closed bottom quark loop. The dashed 
orange lines indicate cuts contributing to the b — > s'yg process. Four-particles cuts of these 
diagrams are kinematically not allowed. 

we use to regulate the ultraviolet, infrared and coUinear singularities. We also stress that 
we used the Feynman gauge for the gluon propagator. 

Since we are interested in the photon energy spectrum, we fix the photon energy by 
inserting in the phase-space integrand on the r.h.s. of (3.1) a factor 

d(^E^-Pi^^ =2mbS{{pb + p.y)''-{l + z)ml) . (3.2) 

Then, after performing the phase-space integrations over the full range, Fi.^^ coincides with 
dG\:^"'{z, fi)/dz given in (2.4) up to a constant factor. 

The optical theorem relates the b quark decay rate we are interested in to the imaginary 
parts of three-loop b quark self-energy diagrams. The contribution of a specific physical 
cut of the three-loop diagrams to the imaginary part of the b quark self-energy can be 
evaluated using the Cutkosky rules [28] , resulting in precisely the same integrals that appear 
on the r.h.s. of (3.1). The various contributions to dG]'^"' {z , fi) / dz can therefore be easily 
represented in terms of cut three-loop diagrams. Figs. 2-5 show the physical cuts that give 
non- vanishing contributions to dG],:f^^{z, fi)/dz. 

The individual contributions to the decay rate (or, equivalently, each cut-diagram) can 
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Figure 5: Same as in Fig. 2 for abelian diagrams without a closed quark loop. The color 
code is the same as in the figures 2, 3 and 4. 
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now be evaluated with the technical tools usually employed in the calculation of multiloop 
Feynman diagrams. To this end, we rewrite all delta functions appearing in the phase-space 
representation (3.1) as well as the one from the kinematic constraint in (3.2) as a difference 
of propagators [29,30], 

K<l" - = ^ { 2 \ -n - 2 \^-n ) ■ (3.3) 
^ ' 2m \q'' — — lO — + lO J 

The expressions for Fi^^ can then be reduced algebraically by means of the systematic 
Laporta algorithm [31], based on the integration- by-part identities (IBPs) first proposed 
in [32,33]. Basically, the Laporta algorithm allows one to express the various contributions 
to the integral in (3.1) (i.e. the various cut diagrams in Figs. 2-5) as linear combinations 
of a small number of simpler integrals, usually referred to as the master integrals (Mis) of 
the problem. The coefficients in front of the Mis are ratios of polynomials of the kinematic 
variables present in the problem; in the case under study, the only kinematic variable is z. 
The reduction procedure is simplified by the fact that all the integrals in which one of the 
cut propagators disappears from the integrand, or is raised to a negative power, are zero, 
because in this cases the ±iO prescription becomes irrelevant. A further comment concerning 
the application of the IBPs to the calculation of phase-space integrals is in order. It is not 
often stressed that the propagators originating from phase-space integrations via (3.3) are 
always accompanied by a factor 6{qo), where qo is the energy flowing in the propagator. In 
building the IBPs it is necessary to take the derivative of the integrand of a given diagram 
with respect to the integration momenta. When the derivative d/dqo acts on 6'(go), the latter 
transforms into a delta function and hence terms like 

are produced. Obviously, the ±iO prescription is of no relevance here and the two terms 
in the parenthesis cancel. Thus terms where the derivative acts on theta functions do not 
contribute to the IBPs. Delta functions enforcing kinematic conditions, like the one in (3.2), 
are not multiplied by theta functions. We made use of two automatic implementations of 
the Laporta algorithm: AIR [34] based on Maple [35] and an independent Mathematica [36] 
code. 

Once the reduction to Mis has been completed, the propagators introduced via (3.3) are 
reconverted into delta functions and the Mis are evaluated. After introducing appropriate 
parameterizations for the three- and the four-particle phase-space [37, 38] (and Feynman 
parameters in order to perform the loop integration if the MI includes a closed loop without 
cut propagators), the simpler integrals can be calculated analytically after expanding in e. 
To obtain analytic expressions for the more difficult ones, it is convenient to employ the 
differential equation method [39] as described in detail in [29,30]. As an independent check 
we also computed all the Mis numerically by means of the sector decomposition method [40] 
(see [9] for further details). 

The color factors for the individual diagrams shown in Figs. 2-5 are listed in Table I of [8] . 
The contribution of the diagrams which are not left-right symmetric must be multiplied by 
a factor 2. Notice that we also included diagrams which arise by replacing the quark loop in 
Fig. 3 by a gluon loop. The four-particle cuts of these diagrams can be calculated in two ways: 
(a) using the complete expression for the polarization sum of the two transverse gluons, or. 
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(b) replacing the polarization sums of the gluons by the negative of the metric tensor; in this 
case, the emission of ghost particles (corresponding to four-particle cuts through diagrams 
with a closed ghost loop) has to be taken into account. We convinced ourselves in [9] that 
the two procedures lead to the same result. In the present paper we used procedure (b). As 
pointed out in [8], the contributions of the diagrams with a closed gluon/ghost loop can be 
obtained from those involving a closed light quark loop by replacing (in the Feynman gauge) 

TnN,^-C^(^l + '- + 0{e')^ . (3.5) 

The diagrams in Figs. 2-5 involve IR, coUinear and UV singularities, that appear in the 
intermediate steps of the calculation as poles in the dimensional regulator e. Some of the 
cuts generate not only single, but also double poles; the latter, however, cancel in the sum 
of all the cuts belonging to a given b quark sclf-cncrgy diagram. The sum of all diagrams 
has a single pole that gets removed when adding the counterterm contributions. 

We conclude this section by summarizing the renormalization procedure adopted; it 
coincides with the one employed in [8,9]. As explained after (2.6), our calculation can be 
restricted to 2; 7^ 1. In this case all counterterms of 0{al) (i.e. those contributing to if^^^ 
defined in (2.10)) are induced by renormalizing the various parameters and wave functions in 
the process b s'-fg. To get the counterterms proportional to the color factors Ca and T^N^, 
only the renormalization of the strong coupling constant is needed. For the term proportional 
to TrNh we need the renormalization of the strong coupling constant and that of the gluon 
wave function. In order to renormalize the contribution proportional to Cp, one needs also 
to consider the renormalization of the operator O7, the renormalization of the bottom and 
strange quark wave functions, and of the bottom quark mass. The renormalization constants, 
which are only needed to 0{as) precision in our calculation, are fixed as follows: 



i) the strong coupling constant is renormalized using the MS scheme: 

= l+(T«A^^i^W + CaK(2)) ^+Oial) , = ^ , K^^^ = ; (3.6) 



ii) the b quark mass contained in the operator O7 and the Wilson coefficient Cf^ are 
renormalized in the MS scheme in our calculation. In the present application we only 
need the product of these two renormalization constants which is given by 

Z^zW-l + CpK,r^ + 0{al), i^rr = ^ ; (3.7) 



iii) all the remaining external fields and the b mass are renormalized in the on-shell scheme. 
The on-shell renormalization constants for the b quark mass and wave function are given 
by 

Zg^^l + CFK,,^ + 0{a',), K2,^Krn,. (3.8) 
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For the wave-function renormalization factor of the gluon-field we find 
Z^^^l + TnNHK,^ + 0{al), K,^-^-\l,. (3.9) 

TT 0£ O 

The s-quark field renormahzation constant is equal to one. [Notice that in our approach 
the effects of the diagrams in Fig. 4 are taken into account through the renormalization 
of the gluon wave function.] 

Consequently, we find {z ^ 1) 

^(2,a) _ ^(2,a,ba.e) ^ ^27^77 + K2b) H^^''^ - Km.H^^'"'''^ , 
^(2,na) ^ ^(2,na,bare) _^ K^) h'^^'^\z) , 

^(2,Nl) ^ ^(2,NL,bare) _^ K^'^ H^^'^\z) , 

^(2,NH) ^ ^J^il) ^ ^^^) ^ ^3_^Q) 

where the quantities H^'^'^'^'"""'' [i = a, na, NL) are the unrenormalized contributions to the 
photon energy spectrum corresponding to the diagrams given in Figs. 2, 3 and 5. The 
function ^(2,NH,baro) yaj^ishes in our approach as just explained. The functions H^^'^^ and 
}jO-,m,e) lY^chide terms of order e and their explicit expressions are collected in the Appendix. 



4 Summary 

A precise calculation of the photon energy spectrum in the inclusive B — > decay is crucial 
for comparing theoretical predictions with measurements. In this paper we calculated the 
NNLO QCD corrections to the photon energy spectrum induced by the magnetic penguin 
operator Of. This subset of contributions is likely to provide the leading correction to the 
photon energy spectrum at 0{al). The present calculation represents an independent check 
of the results of [10], with which we find complete agreement. 

The numerical relevance of the results and some application have already been studied 
in [10] for the normalized photon energy spectrum: in particular, the effect of the non-BLM 
corrections is minor with respect to the dominant BLM component [19]. This statement 
also holds for the unnormalized photon energy spectrum dGrj/dz shown in the left frame of 
Fig. 1, as well as for the integrated quantity [dG77 / dz]dz shown in the right frame of this 
figure, provided the decay width is written to be proportional to m6(/x)^m^ as in (2.1). 

Furthermore, we investigated the impact of the change of the bottom quark mass from the 
pole to the kinetic scheme on observables like the first moment of the spectrum (E^) E-i>Eo 
and the fraction of events with photon energy above Eq, R[Eq). Apart from the change 
of the numerical values of the h quark mass, these observables receive the additional shifts 
presented in (2.20) and (2.25). We find that the results in the two schemes are compatible, 
although the pole scheme is expected to have larger higher order corrections. 

In the calculation we employed a number of techniques usually applied to the calculation 
of multiloop Feynman diagrams, such as the Laporta algorithm, the differential equation 
method for the analytic evaluation of master integrals, and sector decomposition as an 
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independent numerical check of the analytically evaluated master integrals. Several technical 
aspects of the calculation were discussed in detail in Sec. 3. The same techniques can now 
be applied to the calculation of the subleading 0(«f) contributions (O7, Og) and (Og, Og) to 
the photon energy spectrum of the inclusive radiative B decay. 
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A Functions iJ^^'^) and 

Here we collect the functions H^^'^^ and //(i'™'*^) including 0{s) terms as required by rcnor- 
malization. For z ^ 1, that is with terms proportional to S{1 — z) dropped, they are given 

by 



In(l-z) 7 1 1 + z 



1-z Al-z 



ln(l -z) + 



7 + z-2z'^ 



(1 + 4£L^) 



l-z V 3 2 



71n(l-^) 3\n^{l-z) 
4 l-z ^2 l-z 



+ 21^ + 4^^ + ^) ^^^^ 4 ^^^^^ -z) + 2\nz) 

Hl^z)(^-U.il-z))+l±^ 



(Al) 



H^^'"'''\z) 



ln(l - z) 



6 



+ (1 + z) ln(l - - (6 + 2z) 



+e 



l-z l-z 

2 9\t?{1-z) _ h\{l-z) 



il-zY 2 l-z 



l-z 



- 12- 



(1 + 4£L^) 
In z 



l-z 



l-z 



+ (1 + 2) 



TT 



-Uoil-z] 



l-z 2 



- - il + z)hiHl - z) 



+ (6 + 2) ln(l -z) + (12 + Az)\nz-Q-2z 



(A.2) 
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